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FOREWORD 



This review was written to fulfill the requirements for the "Habilitation" pro- 
cedure at the University of Bonn. It summarizes a part of the research I have 
conducted over the last years in the field of microquasars and it should stimulate 
students' interest in this challenging and fascinating topic. 

The astronomical methods used in the range of gamma-rays, X-rays, optical 
and radio wavelengths are first reviewed and then, all of them, directly applied 
always to the same astronomical source, the periodic microquasar LS I +61°303. 

Reviewing the literature in that field, this presentation is of course biased to- 
wards my own contributions for the reason mentioned in the beginning. There- 
fore, I want to acknowledge here all colleagues working like me for years on 
LS I +61°303, first of all the two discoverers of the source: P.C. Gregory and A.R. 
Taylor, but also D. Crampton and J.B. Hutching for the optical observations, F.A. 
Harrison, D.A. Leahy and M. Tavani for the high energy, K,M.V Apparao, and 
R.K, Zamanov for multiband studies, and of course all my collaborators: J.M. 
Paredes, M. Ribd, J. Marti, S. Garrington and M. Peracaula. 



Maria Mas si, December 2004 
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1 INTRODUCTION 



Since the beginning of the 1980s objects like radio-galaxies, quasi-stellar radio 
sources (Quasars), Seyfert Galaxies are simply classified as Active Galactic Nu- 
clei (AGN), because the "energy-engine" is thought to be the same for all of them: 
A super-massive black hole of millions of solar masses accreting from its host 
galaxy (Fig.[T]). AGN with radio-emitting lobes or jets are called radio-loud, the 
others radio-quiet (Ulrich et al. 1997). 

The class of X-ray binary systems is very similar, the "energy-engine" is a 
compact object of only a few solar masses accreting from the companion star 
(Fig. [T]). Up to now there are known 280 X-ray binary systems (Liu et al. 2000, 
2001), but only 18 of them (Fig. EJ) show evidence of a radio-jet and therefore are 
radio-loud applying the same definition as for the AGN. 

The radio-loud subclass of X-ray binary systems includes together with the 
microquasars - objects where high resolution radio interferometric techniques like 
VLBI have given direct evidence of the presence of collimated and relativistic jets 
(Mirabel et al. 1992) - also unresolved radio sources with a flat spectrum, which 
give indirect evidence for continuous ejection. 

I here review how the three important basic components of a microquasar (Fig. 
13) - a compact object, an accretion disk and a collimated relativistic jet - have been 
observed in gamma-rays, X-rays, optical and radio emission. After a basic intro- 
duction of the accretion-ejection processes presented in section |2l the following 
sections describe the astronomical methods: First, presenting their theory and, 
afterwards their application on the source LS I +61°303. In detail: Optical ob- 
servations, reported in Sect. |3j reveal the nature of the compact object: Neutron 
star or black hole. X-ray observations, discussed in Sect. 4, deliver information 
on the accretion disk, while radio observations allow to study the jet; the section 
describes how the simultaneous use of X-ray and radio tools allow to study "the 
disk-jet" connection. The application to LS I +61°303 shows the limits of an 
approximation of the accretion theory which assumes a constant velocity of the 
accretor along the orbit. Section 5 shows that the observational results can be ex- 
plained if one takes into account the geometry of the orbit. The better understand- 
ing of the physical processes motivates for new observations at higher energy. 
Gamma-ray observations are discussed in Sect. 6. The characteristics of the jet 
(morphology, velocity, etc) are derived from high resolution radio-astronomical 
observations and their typical procedures are described in Sect. 7. Finally, the 
conclusion of this multiband approach applied on LS I +61°303 are presented in 
Sect. 8 
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The Class : 

AGN 

energy-engine: 
supermassive BH 



Quasars 




Figure 1 : : The AGN (Top) and X-ray binary classes (Bottom). The energy 
engine for the AGN is a super-massive black hole (> 10 6 M©) accreting from its 
host galaxy. The X-ray binaries are stellar systems formed by a normal star and a 
degenerate object (a neutron star or a black hole of a few solar masses) accreting 
from the companion star. If there is evidence for a radio jet the X-ray binary 
system is defined as radio-loud. 
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Where are the other objects like SS433 ? 

Margon 1980,1984 



Radio-loud X-ray Binaries: 

MICROOUASARS: 



SS433 

1E1740-2942 
GRS 1758-258 
Cygnus X-3 



Circinus X-l 
LSI 61303 
GRS1915+105 
GROJ 1655-40 
XTEJ 1748-288 
CI Cam 
LS 5039 
V461 Sgr 
Cygnus X- 1 
Sco X-l 
XTEJ1550-564 
XTEJ 1859+226 



Margon 1979, Spencer 1979 
Mirabel et al. 1992 
Rodriguez et al. 1992 
elongation Geldzahler et al. 1983 
Spencer et al 1986 
Schalinski et al. 1990, 1995 
Stewart et al. 1993 
Massi et al. 1993 
Mirabel & Rodriguez 1994 
Tingay et al. 1995; Hjellming & Rupen 1995 
Hjllming et al. 1998 
Mioduszewski et al. 1998 
Paredes et al. 2000 
Hjllming et al. 2000 
Stirling et al. 2001 
Fomalont et al. 2001 
Hannikainen et al. 2001 
Brocksopp et al. 2002 



INVERTED /FLAT SPECTRUM COMPACT SOURCES: 

GX339-4 Fender et al. 1997, 

recently resolved by Gallo et al. 2004 
XTEJ1 1 18+480 Fender et al. 2001 



Figure 2: : Radio-loud X-ray binary systems. The radio-loud subclass of X-ray 
binary systems includes together with the microquasars -objects where high res- 
olution radio interferometric techniques like VLBI have given direct evidence of 
the presence of collimated jets - also unresolved radio sources with a flat spec- 
trum that gives indirect evidence for continuous ejection. At the top: The his- 
torical sentence of Margon, who discovered SS433, the first galactic miniature of 
a quasar, unique for several years. Traces of elongation were evident in Cygnus 
X-3 already in 1983. 1E1740-2942 was the first object with radio jets to be called 
a "microquasar". The two sources in red LS I +61°303 and LS 5039 (both dis- 
covered by our group) are the only two microquasars coincident with high-energy 
Gamma-ray sources. 



6 




Figure 3: : The basic components of a microquasar: A spinning compact 
object, an accretion disk and a collimated relativistic jet. The compact object here 
is of a few solar masses, while in the extragalactic analog (AGN) the compact 
object is a black hole of millions of solar masses accreting from its host galaxy. 
The compact object in a microquasar accretes from a normal star in orbital motion 
around it. The study of the periodic velocity- shift of the optical spectral lines of 
the companion star allows to determine the mass of the compact object and to 
establish whether it is a neutron star or a black hole. The inner part of the disk 
emits X-rays. The inner radius is three times the Schwarzschild radius, the outer 
radius a factor of 10 3 larger (the figure is not in scale). Due to magneto-rotational 
instabilities a part of the disk is propelled into a relativistic jet, studied at high 
resolution with radio interferometric techiques. In some microquasars, like SS433 
and LS I +61°303 , the jets are precessing. If the precession brings the jet pointing 
towards the Earth the large variable Doppler boosting mimics the variability of 
Blazars and in this case the Microquasar is called Microblazar (Massi 2003). 
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2 THE ACCRETION-EJECTION PROCESS 
2.1 Accretion 

X-ray binaries are stellar systems formed by two stars of a very different nature: 
A normal star (acting as a mass donor) and a compact object (the accretor) that 
can either be a neutron star or a black hole (White et al. 1996). 

Several mechanisms have been proposed to explain the presence of a compact 
object in a binary system, and they principally depend on the mass of the com- 
panion. If the companion is a low mass star (Low Mass X-ray Binary, LMXB) 
the theory assumes, first, the formation of the neutron star/BH, which later cou- 
pled with its companion in a (tidal) capturing process. As a matter of fact several 
LMXB are close to the core of globular clusters or near the center of the Galactic 
bulge (Verbunt & van den Heuvel 1996). The High Mass X-ray Binary (HMXB) 
systems, where the companion has a mass above 5M0, have a galactic disc dis- 
tribution characteristic that of young stars (population I). It is assumed that large 
scale mass transfer has occurred in the system before the supernova explosion: 
When the progenitor star of the compact object had evolved to a Red Giant and 
had filled its Roche lobe, the smaller companion accreted from it to a level that 
it survived after the explosion. In most HMXRB systems the massive companion 
of the compact object is a rapid rotating Be star, whose formation is explained by 
the large amount of angular momentum received together with matter from the 
initially more massive and therefore faster evolving companion (Verbunt & van 
den Heuvel 1996). 

Since the binary pair is in orbital motion around the common center of gravity, 
the matter leaving the companion star has some angular momentum (J), which 
prevents it from directly falling into the accretor. The stream of matter orbits 
the compact object with a radius determined by / and the mass of the compact 
object (Mx). The angular momentum is redistributed by the viscosity: Some of the 
material takes angular momentum and spreads outwards, whereas other material 
spirals inwards. In this way a disk is created from the initial ring of matter (King 
1996; Longair 1994 p. 135). Gradually the matter drifts inwards until it reaches 
the last stable orbit, called "the inner radius" of the accretion disk (/?,„), which for 
a non rotating black hole is approximately three times the Schwarzschild radius 
(r s ): 

2GM X 



c 2 



(1) 



*in = 3 r s ~ 9 (^|) km. (2) 

The viscosity has two effects: Besides the transport of angular momentum it also 
acts like a frictional force resulting in the dissipation of heat. The amount of fric- 
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tion depends on how fast the gas orbits around the compact object; the temperature 
(7] n ) reaches its maximum at the inner disk where it rises up to (Longair 1994 p. 
141) 

r m ^2 10 7 (^)- 1 /4 K . (3) 

Mq 

On the basis of this equation we see that for a microquasar of about 1 solar mass 
the matter around the last stable orbit is heated up to tens of million degrees 
therefore emitting predominantly in the X-ray band. This led to the name 'X- 
ray binaries' for this class of objects. The Earth's atmosphere is opaque at these 
wavelengths; therefore it is understandable that there was an impasse (until re- 
cent developments in X-ray astronomy) in discovering such stellar sources and 
(besides SS433 discovered by chance) their subclass with relativistic jets. On the 
contrary, the temperature of the last stable orbit around a super massive black hole 
of an AGN with a mass of 10 9 MQ is 7} M = \0 5 K. Therefore the emission is in the 
ultraviolet band causing the "blue bump" associated with the "visible" Quasars. 
This last fact is an example that the same laws are applied to both, AGN and X- 
ray binaries, deriving parameters only scaled with the mass (see Heinz & Sunyaev 
2003; Merloni et al. 2003; Falcke et al. 2003). 
The accretion luminosity can be written as: 

L = r\mc 2 , (4) 

where T], the efficency of energy conversion, expresses here how compact an ob- 
ject (with radius R) is: T] = 1 /2{r s /R) (Longair 1994 p. 134). Whereas for a white 
dwarf T] is only 0.0001, for neutron stars T] is 0.1. As a comparison, the release of 
nuclear binding energy occuring in the conversion of four protons into helium has 
an ^ = — 4^ — " = 7*10 . Thus, accretion in neutron stars already is an order 
of magnitude more efficient as an energy source as compared with nuclear energy 
generation (Longair 1994 p. 134). 

However, there is a limit of energy that is possible to extract by accretion: If 
the force generated in the accretion disk by radiation pressure exceeds the gravita- 
tional force of the compact object a further accretion of gas ceases. The expression 
for that luminosity limit, the Eddington luminosity, is (Longair 1994 p. 137; Frank 
et al. 2002) 

L E = 1.3 10 38 (^)ergs- 1 (5) 

Mq 

2.2 Magnetohydrodynamic Jet Production 

In the case of a small vertical magnetic field threading the disk the plasma pressure 
dominates the magnetic field pressure and the differentially rotating disk bends the 
magnetic field lines, which are passively wound up (Fig. |4|) (Meier et al. 2001). 
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Figure 4: : A differential rotating disk drags the field lines of a vertical mag- 
netic field.Meier, Koide & Uchida 2001. 

www.batse.msfc.nasa.gov/colloquia j abstract sjspring03 / presentations j meier.pdf 

Due to the compression of the magnetic field lines the magnetic pressure may 
become larger than the gas pressure at the surface of the accretion disk, where 
the density is lower. At this point the gas starts to follow the twisted magnetic 
field lines, creating two spinning flows. This extracts angular momentum (mag- 
netic braking) from the surface of the disk and enhances the radial accretion. The 
avalanching material further pulls the deformed magnetic field with it and after- 
wards magnetic reconnection may happen (Matsumoto et al 1996). The flux tubes 
open up and reconnect as is known from stellar flares (Massi et al. 2002). 

The thickness of the disc is a fundamental parameter in this magneto-rotational 
process, or better the extent of the poloidal magnetic field frozen in the disc (Meier 
2001; Meier et al. 2001; Maccarone 2004). No radio jet is associated with X-ray 
binaries in High/Soft states (Sect. 4. 1), where the X-ray spectrum is dominated by 
a geometrically thin (optically thick) accretion disc (Shakura & Sunyaev 1973). In 
the contrary numerical results show a jet being launched from the inner geometri- 
cally thick portion of the accretion disc that is present (ADAF/Corona) when the 
X-ray binaries are in their Low/Hard state (Sect. 4.2) (Meyer et al. 2000; Meier 
2001). 

In conclusion, a better understanding of the transition from radio-quiet to 
radio-loud therefore seems to be possible through a better understanding of the 
X-ray states and their switch mechanism. 
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2.3 Strong Magnetic Fieds: the X-Ray Pulsars 

A radio-loud X-ray binary system may contain either a black hole or a neutron 
star with a low (B < 10 10 Gauss) magnetic field. Accreting neutron stars with a 
low magnetic field can give rise to jet production because of the following reason: 
As described in the previous section only for a low magnetic field can the plasma 
pressure dominate and bend the field. On the contrary, jet formation is prevented 
in presence of a strong magnetic field. In the case of B> 10 12 Gauss, the plasma is 
forced to move along the magnetic field lines, converges onto the magnetic poles 
of the neutron star and there releases its energy creating two X-ray emitting caps 
that, in case of a misalignment of the rotation and the magnetic axis, produce X- 
ray pulses (Fig. 13). X-ray pulsars are not associated with microquasars. The lack 
of detected radio emission from X-ray pulsar systems is discussed in Fender et al. 
(1996). 

X-ray puisair 




Figure 5: : Sketch of an X-ray pulsar. 

http : / / Iheawww.gsfc .nasa.gov / users / white / xrb / xray_pulsar.gif 
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3 OPTICAL OBSERVATIONS 



3.1 The Nature of the Compact Object 

The most reliable method to determine the nature of the compact object is the 
study of the Doppler shift of absorption lines in the spectrum of its companion. 
The study of the changing radial velocity during the orbital motion is a technique 
that has been applied for more than one hundred years to measure the masses 
of stars in binary-systems. The same method is applied for systems like X-ray 
binaries, where one component is "invisible". In this case the variations of the 
radial velocity of the normal companion during its orbit are studied. 



L 1 ' 1 ' i ' ' 1 L * ' 1 ' 1 I 




■ ■ * ■ 1 ■ ■ ■ ■ 1 ■ « ■ ■ 1 ■ - ■ ■■ 
0.6 1 1.5 £ 

Relative Phase 



Figure 6: : Amplitude of the radial velocity variations versus orbital 

phase(Filippenko et al. 1999; GRS 1009-45). Using the Doppler shift of spectral 
lines from the companion star orbiting around the compact object, one determines 
the mass function /, lower limit to the mass Mx of the compact object. Notice 
that for a better display, the orbital phase interval 0-1 is repeated twice. 

The amplitude ( ^f c ) of the radial velocity variations (Fig. |6j) of the mass donor 
and the period (P or b) of the system applying Newton's/Kepler's third law define a 
quantity called the "mass function" (Charles & Wagner 1996), which is equal to: 

= P orb Kj = M x sin 3 i 
J ~ 2nG ~ (M x +M) 2 
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where Mx and M are the masses of the compact object and of the companion, 
respectively, i is the angle between the axis of the orbit and the line of sight and G 
is the gravitational constant. 

The mass function alone already provides a lower limit for Mx correspond- 
ing to a zero-mass companion (M=0) viewed at the maximum inclination angle 
(i=90°). In the cases where the inclination i and the mass of the companion M are 
known one can solve for the mass Mx of the invisible object. 

Rhoades & Ruffini (1974), taking the most extreme equation of state that pro- 
duces the maximum critical mass of a neutron star, established the very upper limit 
of 3.2 MO for a neutron star. This absolute maximum mass provides a decisive 
constraint to observationally distinguish between neutron stars and black holes. 

In Fig. |7] a list of some X-ray binaries is given for which both mass function 
and M x is available. All sources below GRS 1009-45, with f(M) = 3. 17M0, can 
be defined black hole candidates on the basis of the mass function alone. On the 
contrary, for cases where f(M) < 3 the determination of inclination and mass of 
the companion is mandatory to determine the type of object. 

It is worth mentioning here that the accumulation of accreted material on the 
surface of a neutron star triggers thermonuclear bursts (see typical profile in Fig. 
[Hi. These are called bursts of Type I. No Type I burst has ever been observed 
from a compact object where optical observations resulted in a mass above 3 M0. 
That fact might confirm that in black holes there is no surface where material can 
accumulate (Narayan & Heyl 2002). In conclusion: Observations of Type I bursts 
give a direct evidence for the existence of a neutron star. 

3.2 The Nature of the Companion Star 

The classification of the X-ray binaries into Low Mass X-ray Binary and High 
Mass X-ray Binary systems leaves unspecified the nature of the accreting object 
and is based on the mass of the companion star (van Paradijs & McClintock 1996). 

A LMXB contains a late type (K,M) low mass donor star. The mass trans- 
fer takes place via Roche lobe overflow (Frank et al. 2002): Material streams 
through the inner Lagrangian point and will orbit the compact object at the radius 
determined by its specific angular momentum. 

In HMXB systems, the companion is an OB star. OB stars have a substantial 
stellar wind (mass loss rates 10~ 10 — lO _5 M0yr _1 ) with a velocity of v Wi „^ ~ 

v escape = J ^jf- — 10 3 km/s. However, matter leaves the star in all directions, not 
only towards the accretor as in the case of Roche lobe overflow. This accretion 
therefore is less efficient (King 1996). The expression for the accretion rate is 
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(Bondi 1952): 

M=15E=4^«! (6) 

rel 

where p w ind is the density, and v re i depends on the velocity along the orbit v or b 
and on the wind velocity (v^). The accretion therefore becomes more efficient 
for denser and slower winds present in Be stars. In these rapidly spinning stars 
together with a high velocity (1000 km s _1 ) low density wind at high latitudes 
there also exists a dense and slow (100 km s _1 ) disk-like wind around the equator 
having a power law density distribution (Waters et al. 1988). 



3.3 LS I +61°303: The Be-Star 

The ultraviolet spectroscopy of LS I +61°303 by Hutchings & Crampton (1981) 
indicates that the primary star is amain sequence B0-B0. 5 star (L~ 10 38 erg sec -1 , 
T e ff ~ 2.6 10 4 K). Its distance is 2.0±0.2 Kpc (Frail and Hjellming 1991). The 
optical spectrum is that of a rapidly rotating star with Vsin i = 360 ± 25 km s _1 . 
The critical rotational velocity for a normal B0 V star is ~600 km s _1 and Be stars 
may rotate at a velocity that does not generally exceed 0.9 of this, i.e. 540 km s _1 
(Hutchings et al. 1979). Therefore the lower limit for the inclination of the orbit 
compatible with these data is 38°. However, Hutchings & Crampton observed 
shell absorptions in the strong Balmer and He I lines. For a disk sufficiently flat 
this corresponds to a large inclination angle (i ~ 90°) (Kogure 1969). The result 
is a range of 38°-90° for the inclination of the orbit for LS I +61°303 . 

LS I +61°303 is the only X-ray binary system showing variations compatible 
with the orbital period at X-rays (Paredes et al. 1997; Leahy 2001), at Gamma- 
rays (Massi 2004; Massi et al. 2004b), at optical wavelengths in both continuum 
(Maraschi & Treves 1981; Paredes & Figueras 1986; Mendelson & Mazeh 1989) 
and lines, (Zamanov & Marti 2000; Liu et al. 2000; Apparao 2000). The most 
accurate value for the orbital period is however from radio astronomical measure- 
ments resulting in 26.4960 ± 0.0028 days (Gregory & Taylor 1978; Taylor & 
Gregory 1982; Gregory 2002). 

The range of the mass for a B-star is 5-18 M©. Fits performed on near infrared 
data by Marti and Paredes (1995) result in an eccentricity of e~ 0.7-0.8 and mass 
in the range 10-18 M0. 

Finally, because of the high eccentricity, an important parameter of the system 
is the phase at the periastron passage, already determined by Hutchings & Cramp- 
ton (1981) and very recently confirmed by Casares et al. (2004) to be (f)=0.2. The 
zero phase by convention refers to the date to=JD 2443366.775, the date of the 
first radio detection of the system (Gregory & Taylor 1978). 

Figured shows a sketch of the system with the compact object travelling (and 
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accreting) through the dense, variable and structured wind of the Be- star along the 
quite eccentric orbit. 

3.4 LS I +61°303: Really a Neutron Star ? 

The composite and variable nature of the spectral features of LS I +61°303, its 
long period and its high eccentricity make it difficult to derive meaningful radial 
velocities for this source. The observations of K c made by Hutchings & Crampton 
(1981) imply a mass function in the range 0.0028< / <0.043 (see also Punsly 
1999). Following the discussion in section l3~T1 we see that the low value of the 
mass function is not enough to establish the real nature of the compact object. 
Knowledge of the values for inclination and mass of the companion star also is 
necessary for that. On the other hand, we have seen in the previous section that the 
inclination has a quite large range (38°-90°) and the possible range for the mass 
(M) of the companion can be 10-18 M©. 

Hutchings & Crampton (1981) have assumed /=0.02, an inclination of about 
70° (sin 3 i=0.8) and M=10 M© and derived M*=1.2 M©. For this reason in the 
literature it has generally been assumed that the compact object in LS I +61°303 
is a neutron star. Only Punsly (1999) discussed the possibility that the compact 
object in LS I +61°303 could be a black hole and he presented a model for the 
high-energy emission based on it. 

As discussed above and also in Massi (2004) the uncertainties in the param- 
eters derived by optical measurements are rather large; changing the inclination 
to i=38° we already get M x =2.5 M0. If we assume a value M=18 M© we ob- 
tain Mx=3.4 MQ. Recently Casares et al (2004) have determined an upper limit 
of /=0.027 which would correpond to M x =3.8 MQ (Massi et al. 2004b). In 
conclusion: Accounting for the uncertainties in inclination, the mass of the com- 
panion and the mass function it cannot be ruled out that the compact object in 
LS I +61°303 is a black hole. 

4 X-RAY AND RADIO OBSERVATIONS 

An X-ray binary system is called transient, if at least one outburst occurs with 
a flux variation of more than 2-3 orders of magnitude greater than the normal 
flux (McClintock & Reimmllard 2004). This outburst, which may last for days 
to months, is directly related to a variation of the accretion disk and therefore is 
distinct from the outburst of Type I discussed in section 13.11 which lasts a few 
seconds. 

Generally Microquasars have been discovered through high resolution radio 
observations immediately performed during new transients. However, LS 5039 
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was discovered (Paredes et al. 2000) on the basis of a cross-identification in cat- 
alogs (in optical, radio, X- and Gamma-rays) without any X-ray outburst calling 
attention to it. Nowadays one observes changes in the X-ray "states" (Fig. FTOb . 
These states are defined by spectral (see below) and timing characteristics (van 
der Klis 2004): As soon as an X-ray binary is discovered in a "Low/Hard state", 
it immediately is observed at radio wavelengths with high-resolution techniques. 

4.1 High/Soft State and Multicolor Disk 

X-ray binaries with a neutron star as compact object may have spectra that are 
completely different depending on whether the magnetic field of the neutron star 
is strong or weak. 

The form of the spectrum of a binary X-ray pulsar, with a surface magnetic 
field > 10 12 G is a fiat hard power-law function with a sharp cut-off above a few 
tens of keV (Tanaka 1997; White et al. 1996). 

The spectrum of an X-ray binary with a weakly magnetized neutron star is 
typically formed by the properties of the accretion disk and the neutron star enve- 
lope. The neutron star envelope contributes to the harder part of the spectrum and 
has a temperature of ~ 2.5 keV. Mitsuda and collaborators (1984), assuming an 
optically thick disk, where the energy generated by viscosity is locally dissipated 
in blackbody radiation, have represented the disk spectrum as a superposition of 
spectra with temperatures varying from a low value T out at the outer edge to a 
maximum Tj n at the inner edge (i.e. at the the inner radius 7?,„) of the disk. This 
is the reason why the disk is generally called a multi-temperature or multi-color 
disk. By means of this model T,„ and Ri„ can be determined through the softer 
part of the observed spectrum (/(£)), which is represented by: 

jU J T ont ^in ^in 

here i is the inclination angle of the disk, D is the distance and B(E,T) is the 
Planck function (Shakura & Sunyaev 1973; Mitsuda et al. 1984; van Paradijs & 
McClintock 1996; Tanaka 1997). 

X-ray binaries known to contain a black hole, proved by measurements of 
the mass function resulting in a mass > 3 M , have spectra with a soft compo- 
nent accompanied by a hard power-law tail (Tanaka 1997). The soft component 
is described by the multicolor blackbody spectrum given above and therefore it 
is associated with the accretion disk around the black hole. This X-ray state is 
defined: High/Soft. In this respect it is quite interesting to compare the different 
values for T ; „ and Ri„ derived in the two cases of neutron stars and black holes 
(Tanaka 1997). In Fig. [TT]the values of R^cos 1 ' 2 / obtained from the fits for ac- 
cretion discs around black holes and neutron stars are collected: The projected 
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inner radius R ; „ of accretion disks around neutron stars always results in values 
of <10km, while the values for black hole binaries all are larger by a factor of 
3-4 than those for neutron stars. This shows that these compact objects are indeed 
more massive than 3M0 as expected following the relationship Rt n M x (Eq. 0). 

The temperature T,„ for disks around black holes is always found to be less 
than ~1 keV, significantly lower than that for disks around neutron stars with 
similar luminosities. Also this difference is understood in terms of the difference 
in the mass Mx of the compact object: T °= ( 1 /M x ) 1//4 keV (Eq. [3l> . 

No blackbody component is present in the X-ray spectra of black hole X-ray 
binaries, which is consistent with the absence of a solid surface in a black hole. 
The second spectral component in black holes in the High/Soft state is a weak 
power-law with spectral index T, defined by the photon flux °cE~ r (Fig. [T2l) . A 
photon index of 2.0-2.5 has been determined by Tanaka (1997) for a sample of 
5 black holes. Esin and collaborators quote a range from 2.2 to 2.7. The recent 
review by McClintock and Remillard (2004) gives a photon index ranging from 

2.1 to 4.8 for 10 black holes. 

The power law component in the X-ray spectra of accreting black holes in 
their High/Soft state is generally interpreted as the result of inverse Compton up- 
scattering of low-energy disc photons by electrons with a power-law or at least 
a hybrid distribution (consisting of both thermal and non-thermal electrons) that 
can be located in coronal regions (possibly flaring) above the disc (Coppi 2000; 
Zdziarski et al. 2001). 

Evidence for the existence of an accretion disc corona comes from systems 
seen almost edge on: The strong central X-ray source (i.e. the inner disc) remains 
hidden behind the disc rim, but X-rays are still seen. The source of emission must 
be quite extended because the eclipse by the companion star is only partial (White 
et al. 1996). The origin of accretion disc coronae is described by buoyancy of 
magnetic fields amplified in the disk (see Miller & Stone 2000 and references 
therein). 

4.2 The Disk- Jet Connection 

Both X-ray binary systems containing weakly-magnetized neutron stars and those 
containing black holes change their spectral shapes (Tanaka 1997). The most 
drastic change is a transition from a spectrum with the accretion disk component 
as that discussed in the previous section (Fig. fT2l . to a spectrum without it and 
showing a single power-law component (Fig. fT3l) . 

In the sample discussed by Tanaka (1997) the photon index (T) in this state 
varies in the range from 1.4 to 1.7 for systems containing black holes and is equal 
to 1.8 for systems with neutron stars. Esin and collaborators (1998) quote the 
range from 1.4 to 1.9. McClintock and Remillard in their recent review (2004) 
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give r=1.5— 1.9 (excluding GRS 1915+105) for 9 black holes. In conclusion, the 
power law in this state is definitely less steep than in the High/Soft state. This 
spectral state, present also in systems with neutron stars, is defined in the literature 
as Low/Hard only for systems with black holes. 

It has been established that when an X-ray binary system is radio-loud and in 
particular with a flat or inverted radio spectrum (i.e. spectral index a > with 
flux density S v a ) then it is always in its Low/Hard state (Fender 2004 and 
references therein). Emissions in the radio band and at hard X-rays are related by: 

Uadio x L°x ( Corbel et al 2003 ; Gall ° et aL 2003 )- 

FigurefT4lshows a multiband monitoring of GX 339-4. At the beginning of the 

observations both radio emission and hard-X ray emission (Hard/Low state) are 

present, whereas the emission in the softer X-ray band is quite weak. When the 

system switches to the High/Soft state, then both radio and hard X-ray emission 

become quenched. 

Finally, as GX 339-4 switches again into its Low/Hard state, radio emission is 
again observed (Fender et al. 1999). From the plot it is clear that the radio and the 
hard X-ray fluxes are strongly anticorrelated with the soft X-rays. 

As stated above the radio emission during a Low/Hard state has a spectrum 
which is flat or inverted. As shown in the sketch of Fig. [15] one can imagine a 
continuous jet in adiabatic expansion (conical jet) as formed by contiguous cylin- 
ders of increasing radii and decreasing B ; . With each cylinder a canonical syn- 
chrotron spectrum is associated with an optically thin part S v^~ p '' 2 (where 
p is the electron energy index) and an optically thick part S V , with the two 
parts of the spectrum joining around Vp ea k ^ 

B (p+2)/(p+4) 

(van der Laan 1966; 

Dulk 1985). Hence, the composite spectrum in case of a prolongated emission 
will have an optically thin part (that of cylinder 1) and a thick part (cylinder 3) 
at the two opposite ends, but will also develop a central part with an intermediate 
or even flat slope (Torricelli et al. 1998; Massi 1999). A flat spectrum therefore 
reveals a continuous jet. That has observationally been proven by direct imaging 
a radio jet in Cyg X-l during its Low/Hard state (Stirling et al. 2001). 

With sufficient sensitivity in the range of a few tenths of keV in the Low/Hard 
state it has been possible to observe a second component due to the accretion disc 
and therefore to measure the inner disk radius. The 7? ; „ resulted to be > 50r s 
(McClintock et al 2001b; McClintock & Remillard 2004). In the High/Soft state 
(sect. 14.11 and Fig. fTTTl . the measured inner radius of the accretion disk corre- 
sponds to Ri n ~ 3r s . Either the inner disk has been removed or during the transi- 
tion High/Soft to Low/Hard the inner part of the disk has made a transition to a 
cooler state, which makes it effectively invisible in X-rays. The transition from 
a soft thermal state to a hard power-law state therefore corresponds to a change 
in the disk structure, evolving from a state mainly characterized by the emission 
from the inner part of the disk to a state, in which this inner-most region has 
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been strongly modified (Tanaka 1997; Belloni et al. 1997). Since the magneto- 
hydrodynamic theory of jet production (Sec. 12.21) assumes a large vertical mag- 
netic field component, the inner-most region is expected now to be geometrically 
thick and therefore strongly different from the geometrically thin (optically thick) 
case of the High/Soft state. However, how in detail is the space within the inner 
radius filled? What is the origin of the power law emission in the Low/Hard state 
and why is it so well correlated to the radio emission? This all is still matter of 
debate (Fig. EH) ( see the reviews by McClintock & Remillard 2003 and Fender 
2004). 

Together with the two extreme states High/Soft (radio quiet) and Low/Hard 
(continuous ejection) there are two intermediate (concerning their hardness) states: 
the Very High State (VHS) and the Intermediate State (IS). Multiple recurrent os- 
cillations in X-rays in the source GRS 1915+105 (Belloni et al. 2000) are due to 
different VHS-like states called A,B and C, where ejections- emitting in the radio 
band - occur in the hardest (C) state (reaching a Low/Hard state with T ~ 1.8). 
Synchronized variations of the inner radius have been observed to occur during 
this oscillations (Belloni et al. 1997,1997b). Figure shows the spectral change 
corresponding to a variation of Rj n from 20 km (value compatible with the last 
stable orbit around a rotating black hole) to more than 300 km. During similar 
episodes of cyclic variations of the inner disk in GRS 1915+105, Mirabel and col- 
laborators (1998) could follow the onset of a flare - first at infrared wavelenghts 
and then at radio wavelengths - with a delay consistent with synchrotron radi- 
ation from expanding magnetized clouds of relativistic particles (FigfT8l. The 
straightforward interpretation was that during the disappearance of the inner disk 
a relativistic plasma cloud was expelled. The mass of the ejected cloud has been 
estimated to ~ 10 19 g (Mirabel 1998), while the matter which disappeared from 
the inner disk in one dip of similar time length (Belloni 1997), has been estimated 
to ~ 10 21 g. This fact could imply that only a very small fraction of the mass is 
ejected, whereas is not clear which fraction fell indeed onto the compact object 
(Mirabel and Rodriguez 1999). 

These small oscillations are now well established, they have recurrence time- 
scales of tens of minutes and have been observed in X-ray (interpreted as possible 
draining and refilling of the inner disk) and at infrared, millimeter and radio wave- 
lengths (interpreted as repeated ejection events) (Fender et al. 1999: Eikenberry 
et al. 1998; Fender et al. 2002). Finally, Fender and coauthors (1999: 2002) 
have shown that these low (~ 40 mJy in radio) amplitude oscillations can happen 
during the decrease of flux in major flare events. 

In conclusion: 

1. A Soft state, characterized by disk emission (at temperatures ~ 1 keV con- 
tributing mostly at 1-2.5 kev) and a power law component steeper than T > 2 
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implies a radio-quiet X-ray binary system. 

2. Quasi-periodic oscillations with time-scales of several minutes in X-ray and 
in the infrared or radio band are a signature of episodic disk-removal/variations 
and plasma bubble ejection. In case of such isolated small ejections, one can 
follow the adiabatic expansion of the cloud and monitor radiation becoming 
optically thin at progressively lower frequencies. 

3. When the X-ray binary system is persistently emitting a radio jet the X- 
ray spectrum has a power law component with T ~ 1.6. This X-ray state 
is called Low/Hard. The superposition (Fig. fT5t of spectra of different 
contiguous jet regions with different self-absorption cutoffs result in a com- 
posite flat spectrum (i.e. S<* v a , with a ~ 0) through and beyound radio 
wavelengths. Emissions in the radio band and at hard X-rays are related by: 

Lradio °^ 

4.3 LS I +61°303: Soft and Hard States 

There are three X-ray observations of LS I +61°303: with ROSAT (Taylor et al. 
1996), ASCA (Leahy et al. 1997) and RXTE (Harrison et al. 2000). The ROSAT 
observation was performed over a total orbital cycle in the energy range from 0.07 
to 2.48 keV. A single component fit was made, either a black body or a power-law 
and average fitted results were presented: A temperature of 0.26 keV and a power 
law index of 1=2. The hardness-ratio is calculated all around the orbit and Taylor 
et al. (1996) notice the hardening of the X-ray emission during the onset of the 
second radio outburst (Fig. fT9l) . 

Greiner and Rau (2001) calculated the photon index around the orbit based on 
RXTE data (measured range 2.3-25 keV) (Fig. l20b . It is quite remarkable that 
the fit (only a single component, i.e. a power law) gives Y =2.0-2.4 all around the 
orbit (i.e. a High/Soft state) except for one point - simultaneous with the onset of a 
radio outburst -, where the photon index is 1 = 1.6 (i.e. Hard state). This seems to 
be the typical case of the theory presented in the previous section for a Low/Hard 
state. The X-ray peak occured at phase ~0.48 whereas the radio outburst came 
almost 6 days later (Fig. l20l . However, in the light of developments in the theory 
of the disk-jet connection (reviewed in the previous sections) the spectral switch 
to the hard state is relevant for the ejection; this switch must precede the onset 
of the radio outburst and RXTE data actually show this for LS I +61° 303 (Mas si 
2004). 

If a source remains in the Low/Hard state, it is radio-loud and the spectrum 
is flat. How does the spectrum of LS I +61°303 behave taking into account that 
the Low/Hard state seems to be of a quite short duration ? Indeed, a flat spectrum 
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has been measured by Taylor et al. (1996): in Fig. [T9Tmiddle the radio spectral 
index is equal to zero before the radio peak. When the ejection phase is terminated 
the figure very nicely shows that the flat spectrum evolves into an optically thin 
one (i.e. a changes from 0, due to the composite spectrum, to -0.3 related to the 
optically thin part of the "last"emitted bubble or cylinder as in Fig. fT5t . As shown 
by Paredes et al. (1991) the flat spectrum in LS I +61°303 can be reproduced by 
a model of an adiabatically expanding cloud of synchrotron-emitting relativistic 
electrons only if a continuos ejection of particles (lasting two days) is taken into 
account as well as adiabatic expansion losses. 

Leahy et al. (1997) reported two ASCA observations, where the photon index 
r was 1.63-1.78 at orbital phase (|)=0.2, which is the periastron passage (where an 
ejection is predicted), and r=1.75-1.90 at orbital phase 0.42 coincident with the 
onset of a radio outburst (as shown in Fig. I2T1) . In conclusions, both ASCA values 
give a photon index consistent with a Low/Hard state during predicted/observed 
ejections. 

Harrison and collaborators (2000) have performed a periodicity analysis of the 
two ASCA observations. The result is a clear periodicity (Fig. l22l) in the ASCA 
pointing related to the onset of a radio outburst while no periodicities are found 
in the other more sparcely sampled pointing. The X-ray oscillations occured at 
phase 0.42 before the radio peak, which (Fig. I2TT) occured at phase 0.5-0.6. The 
photon indexes given above reflect an average of all the data; therefore we cannot 
check, if the Low/Hard state is stable or if there is a continuos switching between 
a sort of Very High States (i.e. A,B,C) like for GRS 1915+105 (Belloni et al. 
2000) or if the Low/Hard state is indeed reached only in the hardest interval. Such 
oscillations are also present in the radio band. Peracaula and collaborators (1997) 
have performed a period-analysis for three radio observations: two in a decreasing 
phase of large outbursts and one at a high, but quiescent flux level. While in the 
last data set there was no evidence for a periodicity, on the contrary, a period of 
P = 84 minutes and significant power also at harmonics of P/2 (i.e. ~ 40 minutes 
again) and 2P have been found in the two data sets related to the decay of radio 
outbursts (Fig.l23t. 

In conclusion: 

1 . There is no estimate of any inner radius for the accretion disk in LS I +6 1 °303 . 
The multicolor disk mostly emits below 2.5 keV and this range is excluded 
by the RXTE analysis. 

2. The available X-ray observations for LS I +61°303 reveal transitions to 
Hard/Low states at the onset of radio outbursts as is expected in the con- 
text of the disk-jet connection (Massi 2004). 

3. A prolongated ejection of particles generates a flat spectrum in the late por- 
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tion of the rise in flux. 



4. X-ray oscillations and radio-oscillations are present at the early rise phase 
and during the decay of radio outbursts (Peracaula et al. 1997: Harrison et 
al. 2000).. 

4.4 The Periodical Radio Outbursts of LS I +61°303 

The greatest peculiarity of LS I +61°303 is its periodic radio outburst activity 
with P=26.496 days (Gregory 2002) (see in Figl24l Left-Top a typical radio light 
curve). The second peculiarity of LS I +61°303 is that the amplitude of each 
outburst is not randomly varying, but itself periodic with a periodicity of 4.6 years 
(Fig. 1241 Left-Bottom) (Gregory 1999, 2002). The orbital phase (j) at which these 
outbursts occur is modulated (Gregory et al. 1999) and varies within the interval 
0.5-0.8 (Fig. ElRight and Fig.[25])(Paredes et al. 1990). 

The theory of accretion (Eq. |6l) predicts maximum accretion where the density 
is highest. The maximum density is obviously at the periastron passage, because 
the density of the wind there is the largest and in addition there probably occurs di- 
rect accretion from the the star (Roche Lobe overflow). As shown in Fig. l24l Right 
the (|) at periastron passage is 0.2 and therefore practically is opposite to the orbital 
region, where the radio outbursts occur. Therefore one of the fundamental ques- 
tions concerning the periodic radio outbursts of LS I +61°303 has been for years: 
Why are the radio outbursts shifted with respect to the periastron passage ? 

5 THEORY OF THE ACCRETION: THE TWO PEAK 
ACCRETION MODEL 

Equation[6]gives only one accretion peak for variable density and constant velocity 
v re [. However, the orbit of LS I +61°303 is quite eccentric and therefore with a 
strong variation of the velocity along the orbit. Taylor et al. (1992) and Marti & 
Paredes (1995) have shown that in this case the accretion rate M develops 

V rel 

indeed two peaks: The first peak corresponds to the periastron passage (highest 
density), while the second peak occurs when the drop in the relative velocity v re i 
compensates the decrease in density (because of the inverse cube dependence) 
(Fig. |26l Top Right). 

This also implies that, while the first peak always occurs at periastron passage, 
the second peak may move to different points in the orbit, if variations in p w i n d or 
v re i occur. Figure[26] Right-Bottom shows how for increasing values of the wind- 
velocity the second peak shifts toward the first peak (Marti & Paredes 1995). 
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Figure 7: : Black hole candidates. Compact objects with a mass (M x ) greater 
than 3 MO, upper limit for a stable neutron star. Ramesh Narayan. 

http://cgpg.gravity.psu.edu/events/conferences/Gravitation_Decennial/ 
Proceedings / Plenaries / Sunday / Narayan / narayan .plenary, pdf 
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Figure 8: :Type I X-ray burst. This type of burst comes from a thermonuclear 
flash on the surface of an accreting neutron star (Tennant et al 1986; Circinus 
X-l). No type I burst has never been observed for accreting black holes. 
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Figure 9: : Sketch of the binary system LS I +61°303 . Orbital phases of the 
compact object travelling through the wind of the Be companion star. The ac- 
cretion rate in an eccentric orbit has two peaks: One peak corresponds to the 
periastron passage (cj>=0.2) because of the highest density; the second peak occurs 
in the phase interval 0.4-0.8 (i.e. around apastron) where the drop in the relative 
velocity of the compact object compensates the decrease in density. For supercrit- 
ical accretion, theory predicts matter be ejected outwards in two jets perpendicular 
to the accretion disk plane. However, near periastron inverse Compton losses are 
severe due to the very near Be star (gamma-ray emission is observed but no radio 
emission). During the second accretion peak, the compact object is much farther 
away from the Be star and both, inverse Compton losses and wind opacity, are 
lower: The electrons can propagate out of the orbital plane and radio outbursts are 
observed. The radio emission has been resolved with VLBI observations and in 
agreement with the theoretical predictions the image shows bipolar jets emerging 
from a central core (Mas si 2004b). 
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Figure 10: The X-ray states of a black hole . Left: The states are defined by the 
spectral and timing properties. Right: Morphology of the accretion disk/Corona 
and radio properties during the X-ray states (Fender 2002). 
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Figure 11: : Best fit of Ri n cos 1 ' 2 i for black holes and neutron stars. From 
Tanaka (1997). Expected values are ~ 10 km and ~ 30 km for a neutron star and 
a non rotating black hole of 3 solar masses respectively. 
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Figure 12: Spectrum of XTE J1720-318 in a Soft state. Cadolle Bel et al. 2004. 
Dashed: Multicolor Disk with R„,=84 km. Dotted dashed: Power law photon 
index 2.6 Thick: Total model. 
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Figure 13: Spectrum of XTE J1720-318 during a Low/Hard state. Cadolle Bel 
et al. 2004. Top: Power law photon index 1.8. Bottom: Residuals 
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Figure 14: Radio (Top), hard (Middle), and soft X-Ray (Bottom) monitoring 
of GX 339-4. Fender et al. (1999). Radio emission and Hard/Low state are 
both strongly anticorrelated with the High/Soft state, where the radio emission is 
suppressed. 
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Figure 15: : Inverted spectrum of an expanding continuous jet. A conical jet 
is here represented by three regions, each with a magnetic field Bj. Each region 
produces a synchrotron spectrum with an optically thin and an optically thick 
part joining at Vpeak 06 Bf 1 (i.e. for electron energy index p=2). The composite 
spectrum from the three spectra develops therefore a central part at intermediate 
slope (Massi 1999). 
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Figure 16: : Sketch of possible disk configurations during the High/Soft and 
Low/Hard states. Fender et al. 1999. In the High/Soft state the inner radius 
reaches the last stable orbit Ri n ~ 3r s and no radio jet is present. In the Low/Hard 
state a radio jet is present and Rj n is larger than before (ADAF/Corona/base of jet, 
are all possible models for what is filling now that space). 
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Figure 17: : The unstable central disk of GRS 1915+105. Top: The 2.0-13.3 
keV PCA light curve and the corresponding hardness ratio (13.3-60 keV/2.0-13.3 
keV) (Belloni et al 1997a). The signal oscillates between a harder (i.e. higher 
hardness ratio) state called here "2" and a softer state ( called "1"). Middle: PCA 
energy spectrum for average 1 (circles) and 2 (points) states (Belloni et al 1997a). 
Table: Best-fit R,„ (km) by Belloni et al. (1997a). Bottom: sketch of the author, 
not in scale. The inner part of the disk has made a transition (1 — ► 2) to a cooler 
state, which makes it effectively invisible in X-rays. Owing to approximations 
in the disk blackbody model the value of the inner radius might be unaccurate 
(Merloni et al 2000), however such a large variation cannot be due to these ap- 
proximations (Fender & Belloni 2004). 



32 




0.2- 

0.18 — 



"V/ 



Seat W, J 99? J 



■J I u 



5 o.ofif 



i i i i I i ' '' r > » 
|i| fl.S p 

UT Time (hours) 



fl.fi 



fl.fl 



Figure 18: : Accretion instabilities and jet formation in GRS 1915+105. 

Mirabel et al. 1998. The infrared flare starts during the recovery from the X-ray 
dip, when an X-ray spike was observed. These observations show the connection 
between the rapid disappearance and follow up replenishment of the inner accre- 
tion disk seen in X-rays, and the ejection of relativistic plasma clouds observed as 
synchrotron emission at infrared and radio wavelengths. The hardness ratio (13- 
60 keV)/(2-13 keV) is shown at the bottom. 
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Figure 19: : X-ray and radio flux density and spectral variations of 
LS I +61°303 .Taylor et al. 1996. Top: The 4.9 GHz radio flux density (solid 
line and circles) and the X-ray flux (dashed line and triangles). Middle: Radio 
spectral index. Bottom:Hardness ratio defined as the ratio of photon counts in the 
energy range 1.0-2.48 keV to those in the 0.07-1.0 keV range. 



34 




Figure 20: : Variations of the X-ray spectral fit parameters (top two panels) 
with Time. Greiner & Rau 2001 . The errors correspond to 3a. The bottom panels 
show the count rate from 2.3-25 keV and the radio flux (Green Bank Interferom- 
eter). The photon index (Top panel) drops to the value of 1.6 at JD (2440000.0+) 
10169, in coincidence with the onset of the radio outburst. 
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Figure 21: : X-ray and radio monitoring of LS I +61°303 . Harrison et al. 
2000. Two-frequency radio light curves (bottom) taken simultaneously with the 
two ASCA pointings (flux levels shown in the top panel). Both ASCA pointings 
have photon indexes 1 = 1.5-1.9. The first pointing is at orbital phase (|) =0.2 
which is the periastron passage, where super-accretion is predicted. The second 
ASCA pointing is at the onset of a radio outburst (bottom panel). 
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Figure 22: : Short-term variability in LS I +61°303 . Harrison et al. 2000. 
Count rate in the 15 keV band for the ASCA pointing simultaneous with the radio 
outburst (see text). 
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Figure 23: : Radio microflares. Peracaula et al. 1997. Top: Radio oscillations 
with mJy amplitude. Bottom: Mean ra djp light curve of LS I +61°303 obtained 
by folding with a period of 1 .4 h. For a oetter display the orbital phase interval 
0-1 is repeated twice. 




Figure 24: : Radio outbursts in LS I +61°303 . Right: Orbital phase interval (in 
yellow) where the periodical radio outbursts occur. Left-Top: Typical radio light 
curve during a strong outburst.Left-Bottom: The 4.6 yr periodic modulation of the 
amplitude of the outbursts. 



39 




teuuuuuuuuiiM 



Figure 25: : Average radio light curves. Paredes et al. 1990. The average curves 
are relative to the maxima (dash-dotted) and to the minima phases of the long- 
term modulation. Low outburst have a broad distribution 0.4-1.0, whereas strong 
outbursts occur in the range of orbital phases 0.5-0.8. . 
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Figure 26: : The accretion model for an eccentric orbit. Marti & Paredes 
1995. Top: Accretion rate versus stellar wind. The vertical axis is in units of the 
Eddington accretion limit, whose limit is indicated by the dashed line. Bottom: 
Accretion rate for different velocities of the stellar wind. The values are in km s _1 . 
Note how the second super-critical event shifts gradually towards earlier orbital 
phases for high values of the wind velocity. 

On the other hand variations in the mass loss of the Be star have been well 
established by Ha emission line observations (Zamanov & Marti 2000). Gregory 
& Neish (2002) suggest a periodic outward moving density enhancement (i.e., 
shell ejection) in the Be star wind. The variation of Ha emission line (Zamanov & 
Marti 2000) is periodic with a comparable scale (1584 d) as the radio modulation 
and it is in phase with the onset of the outbursts (Fig. l27t . The orbital shift 
in the phase of the radio outbursts is therefore related to variations of the wind 
parameters. 

Finally, Marti & Paredes (1995) have shown that during both peaks the ac- 
cretion rate is above the Eddington limit and therefore one expects that matter is 
ejected twice within the 26.496 days interval. 

In conclusion, radio outbursts displaced from periastron passage correspond 
to the second peak of the two-accretion/ejection peaks. The remaining problem 
therefore is: Why is the first outburst at periastron passage in the radio band miss- 
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Figure 27: : Correlation between Ha emission line variations and the radio 
outburst orbital phase variation. Radio and Ha parameters folded on the long 
(here ~ 1584 d) period of modulation of the amplitude of the radio outbursts. Top: 
Total equivalent width of the Ha emission line Bottom: Averaged values of the 
onset phase. 

ing? 
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6 GAMMA-RAY OBSERVATIONS 

6.1 EGRET Sources 

Third EGRET Catalog 

E> 100 Me V 

+9Q 
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+ Active Galactic Nuclei ■ Pulsars 

• Unidentified EGRET Sources LMC 

• Solar FLare 



Figure 28: Third EGRET catalog, http://cossc.gsfc.nasa.gov/egret 



The Third EGRET Catalog contains about 170 not yet identified high energy 
Gamma-ray sources (E > 100 MeV) ((Fig. l28t . The discovery of the coincidence 
of the microquasar LS5039 (Fig. l29l) with an unidentified EGRET source by Pare- 
des and collaborators (2000) has opened the possibility that other EGRET sources 
could be microquasars as well. Gamma-rays can be produced by external Comp- 
ton scattering of stellar UV photons of the massive companion by the relativistic 
electrons of the jet. LS5039 is a persistent radio emitting source and the Gamma 
ray flux, with all uncertainties reflected by the poor sampling, still reflects this 
persistence (Fig. l29l)(Paredes et al 2000). Therefore, for a periodic source like 
LS I +61°303 periodic Gamma-ray emission would be expected. 
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Figure 29: A high-energy gamma-ray-emitting persistent microquasar: Pare- 
des et al 2000. Right: The location EGRET-map of Hartman et al. (1999). The 
contours are not intensity contours, but statistical ones representing the 50%, 68%, 
95% and 99% probability that a single source lies within the given contour. The 
red dot inside the 95% confidence contour, whose radius is about half degree, is 
the position of LS 5039. Left: VLBA-map of LS 5039. The presence of radio 
jets in this high mass X-ray binary is the main evidence of an accretion process 
resulting in the ejection of relativistic particles. The overall length of the source is 
about 18 AU. Bottom: EGRET fluxes and radio observations showing a roughly 
persistent level of emission over a decade in both bands. 
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LS I +61 303 AS A GAMMA-RAY SOURCE 




136.(1 135.5 135.1) 

Galactic Lcngiliide 



Figure 30: Contour plot in galactic coordinates of 2EG J0421+6119. (Kniffen 
et al. 1997) see text. 

6.2 The Variable Gamma-Ray Source LS I +61°303 

Gregory and Taylor (1978) reported the discovery of a radio source (later on as- 
sociated with LS I +61°303 ) within the la error circle of the COS B y-ray source 
2CG 135+01. This association however remained controversial because of the 
presence of the quasar QSO 0241+622 within the relatively large COS B error- 
box. In the CGRO mission (1991 May- 1995 October) the source, given there as 
2EG J0241+6119, was detected by EGRET with a significance of 17 a, with a 
time averaged photon flux of 9.2 ± 0.6 x 10~ 7 cmT 2 s _1 for energies > 100 MeV. 
This flux reported by Kniffen et al. (1997) is slightly different from those quoted 
in the EGRET Catalogs because of the addition of data after the 1993 September 
cutoff date for the catalog. The position of this gamma-ray source is / = 135°. 58, 
b = 1°.13. As shown in Fig.|30|the contour position obtained with additional data 
is about 1 1 ' from the 2EG catalog position and about 40' from the old 2CG cata- 
log position. The radius of the 95% confidence error contour is about 13', ruling 
out the possible identification with QSO 0241+622 at 1 = 135°.7, b = 2° .2, which 
is 64' away. The position is only 8' distant from LS I +61°303 (Kniffen et al. 
1997). In 1998 Tavani and collaborators established the possibility of variability 
of 2CG 135+01 on timescales of days (Tavani et al. 1998). Massi (2004) exam- 
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ined the EGRET data as a function of the orbital phase and noticed the clustering 
of high flux values around periastron passage. Figure|3T]shows (Massi et al. 2004) 
the follow-up of the EGRET gamma-ray emission along one full orbit. At epoch 
JD 2 450334 (i.e. circles in the plot, with empty circles indicating upper limits) 
the orbit has been well sampled at all phases: A clear peak is centered at periastron 
passage 0.2 and 1.2. At a previous epoch (JD 2449045; triangles in the plot) the 
sampling is incomplete, but the data show an increase near periastron passage at 
§ ~0.3, and a peak at (|) ^0.5. The 3 squares refer to a third epoch (JD 2 449 471). 
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Figure 31: : EGRET data vs. orbital Phase Massi et al 2004b (see text). 

In conclusion, the gamma-ray data strongly support the ejection at periastron 
passage predicted by the two-peak model (Massi 2004b). During the first ejection 
(because of the proximity of the Be star) stellar photons are upscattered by the 
inverse Compton process by the relativistic electrons of the jet (Bosch-Ramon 
& Paredes 2004). The inverse Compton losses are so severe that no electrons 
survive: radio outbursts indeed never have been observed at periastron passage in 
more than 20 years of radio flux measurements (Gregory 2002). At the second 
accretion peak the compact object is far enough away from the Be- star, so that 
energetic losses are smaller and electrons can propagate out of the orbital plane. 
At this point the gamma-ray peak at <j> ~ 0.5 is very interesting. It could originate 
from a second ejection which occurred still enough close to the Be-star. In fact, 
while the first ejection is always at periastron passage, the second ejection occurs 
at a varying point in the orbital phase interval 0.4-0.8. 
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7 RADIO INTERFEROMETRY: IMAGING AT HIGH 
RESOLUTION 



Nowadays it is possible to obtain images of jets at infrared wavelengths and in 
X-rays (Sams et al. 1996; Corbel et al. 2002). However, these jets at tenths of 
arcseconds are not related to the emitting regions close to the engine (with quite 
short lifetimes because of their large adiabatic/synchrotron losses) but require a re- 
acceleration mechanism. The study of the jet closest as possible to the "engine" 
at a spatial resolution up to milliarcseconds (mas) is possible at radio wavelengths 
thanks to Very Long Baseline Radio Interferometry (VLBI) (Appendix). 



7.1 The Jet Velocity 

For symmetric ejection of two jets at a velocity (3 (i.e. expressed as a fraction of 
c), the two (approaching and receding) jets move with an apparent velocity of p\.r 
(Mirabel & Rodriguez 1994, Fender 2004): 

_ PsinO 

pa ' r -i T pcose' (8) 

is the angle between the direction of motion of the ejecta and the line of sight to 
the observer. Depending on the angle, for a jet with (3 > 0.7 the apparent velocity 
(3 a of the approaching jet can become greater than 1 (superluminal effect, see 

Fig.m. 

In order to show how the apparent velocity of the jet is derived, let us assume 
0=90. In this case the proper motion, ^, of the jet on the sky plane is: 

" = TT- 

where the distance D is in kpc and fj is expressed in milliarcseconds per day 
(mas/day). The range of |3 is about 0.15-0.99 and the range of the distance is 
about 1-12.5 kpc. From the two extreme conditions, i.e j^il anc [ J^-, the proper 
motion ranges from 2 mas/day to 170 mas/day. In order to estimate [3 from multi- 
epoch observations one has to select the proper radio network. One must take 
into account that at ?i=6cm the VLBI provides a resolution of ~ 1 mas, MERLIN 
one of ~ 50 mas and the VLA in the largest configuration one of ~ 100 mas. 
Therefore high proper motions are best studied with MERLIN and the VLA. 

Beside multi-epoch observations even with only one observation it is possible 
to recover the quantity [3cos0, if is significantly less than 90°. In fact, the 
observed flux densities S a ,r from the approaching and receding jets, 

S a , = S8 k a - a , (10) 
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Figure 32: Superluminal motion: Apparent transverse velocity for the approach- 
ing jet as a function of angle for different values of the true velocity, from (3=0.63 
to (3=0.99 by steps of 0.03. 



(where a is the spectral index of the emission S <x v a and k is 2 for a continuous 
jet and 3 for discrete condensations) are governed by the Doppler factor, 

(id 



^ y(lTpcos0)' 

(where y= (1 — (3 2 )~ 1//2 is the Lorentz factor) and therefore (Mirabel & Rodriguez 
1994) one can determine the quantity [3cos0 by means of the ratio between the 
flux densities from the approaching and receding jet: 

S a _ fl+f3cos0y- a 



S r Vi-pcosO 

Let us assume an ejection nearly aligned to the line of sight with 0=0.5 and 
with the other parameters: a = —0.5, k=2 and |3 ~ 0.6. Using equation ITOl one 



48 



determines 8^ a — 6 and 8^ ~ a = 0.1. As a result, the counter-jet can be rather 
faint, and if S r results to be lower than the noise limit of the radio image the 
counter-jet will completely disappear. In this case the image will show a one- 
sided jet (the approaching one) and only a lower limit for (3 cos can be estimated 
using the noise limit of the image (Massi et al. 2001). 

A constant ejection angle implies a constant ratio between the flux densities 
from the approaching and receding jet during the epochs. An obvious variation 
of this ratio is interpreted as a variation of the ejection angle , explained as jet 
precession. 

7.2 The Processing Jet of LS I +61303 

The first VLBI observation resolving the source (Massi et al. 1993) in Fig. |33] 
reveals a complex morphology (Fig |33] Bottom): A structure at PA~ 30° formed 
by two components separated 0.9 mas (about 2 AU at the distance of 2 kpc) is 
surrounded by an envelope clearly rotated with respect to it. This envelope could 
be an older expanding jet, previously ejected at another angle (because of preces- 
sion). 

Taylor and collaborators (2000) performing VLBI observations in combina- 
tion with the HALCA orbiting antenna of the VSOP mission, imaged a structure 
reminiscent of the precessing radio jet seen in SS433 (Fig. |34|). On the other 
hand at a lower resolution with a scale up to tens of AU (with the EVN), Massi and 
collaborators (2001) obtained an image that for the first time showed an elonga- 
tion in a clear direction without any ambiguity (see Fig. l33t . The most interesting 
aspect of the EVN map is that, for the first time, we detected asymmetric emis- 
sion in the southeast direction. Using the noise level (a) of the map and the peak 
value of the approaching component (5a 6ak ) for k=2 and a=-0.5 we determine 
[3cos0 > 0.6. This would correspond to the two limits of < 53° for |3 ~ 1 and 
P > 0.6 for = 0°. A value for the velocity well within the range 0.1 c to 0.9 c 
found for other Microquasars (Mirabel & Rodriguez 1999i). 

Two observations at still lower resolution have been performed with MERLIN 
(see Table 1). The first MERLIN image shows a double S-shaped jet extending to 
about 200 AU on both sides of a central source. The morphology of the MER- 
LIN image (Fig. l35k) has a bent, S-like structure. In the small box in Fig. l35h we 
show the simulated radio emission from the Hjellming & Johnston 1988 model of 
the precessing jet of SS 433 (rotated here for comparison purposes). The similar- 
ity between the MERLIN image of LS I +61°303 and the precessing model for 
SS 433 strongly suggests a precession of the jet of LS I +61°303. The preces- 
sion becomes evident in the second MERLIN image, shown in Fig. l35b. where a 
new feature oriented to Northeast at a position angle (PA) of 67° is present. The 
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Figure 33: : The relativistic jet of LS I +61°303 Left-Top: EVN uniform 
weighted map of LS I +61°303 at 6 cm. Only the approaching jet is visible, 
whereas the receding jet is attenuated below the sensitivity limit of the image. 
One-sided jets are the observational evidence that the angle between the ejecta 
and the line of sight is smaller than 90° and that the jet is relativistic. Left- 
Bottom: VLBI observation of LS I +61°303 at 6cm. A structure at PA~ 30° 
(and therefore rotated in respect to the one-sided jet of the EVN image) formed 
by two components separated about 2 AU is surrounded by an envelope clearly 
rotated with respect to it. This envelope could be an older expanding jet, previ- 
ously ejected at another angle (because of precession). Right: The stellar system. 
The accretor, whose disk is "feeding" the jet (mapped with EVN and VLBI) and 
the companion Be-star with its equatorial disk. 
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Figure 34: : VLBI observations in combination with the HALCA orbiting 
antenna.Taylor et al. 2000. 



Date 


Start MJD 


Stop MJD 


^start 


^stop 


April 22 


52021.73 


52022.10 


0.670 


0.684 


April 23 


52022.68 


52023.17 


0.706 


0.724 



Table 1: Log of the MERLIN observations. Start and Stop are given in Modified 
Julian Date (MJD=JD— 2400000.5). The corresponding orbital phases have been 
calculated using the new ephemerides, to=JD 2443366.775 and i , =26.4960 d, from 
Gregory 2002). 
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Figure 35: : The precessing jet of LS I +61°303 . Massi et al 2004. a) MERLIN 
self-calibrated image of LS I +61°303 at 5 GHz using natural weights, obtained 
on 2001 April 22. North is up and East is to the left. The synthesized beam has 
a size of 51 x 58 mas, with a PA of 17°. The contour levels are at -3, 3, 4, 5, 
6, 7, 8, 9, 10, 20, 40, 80, and 160a, being a=0.14 mJy beam" 1 . The S-shaped 
morphology strongly recalls the precessing jet of SS 433, whose simulated radio 
emission (Fig. 6b in Hjellming & Johnston 1988), (rotated here for comparison 
purposes) is given in the small box. b) Same as before but for the April 23 run and 
using uniform weights (see text). The synthesized beam has a size of 39 x 49 mas, 
with a PA of 5 —10°. The contour levels are the same as those used in the April 
22 image but up to 320a, with a=0.12 mJy beam -1 . 
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Northwest-Southeast jet of Fig. l35k has PA=124°. Therefore, a quite large ro- 
tation has occurred in only 24 hours. This fast precession causes a deformation 
of the morphology during the second observation, and the one-sided jet appears 
bent in Figl35b. Only 3a features can be associated with the double jet of the day 
before. The feature at 3a to the East is well compatible with a displacement of 
0.6c x 24 hours. 

The appearance of successive ejections of a precessing jet with ballistic mo- 
tion of each ejection is, as shown in Fig. |36l a curved path that, depending on 
the modality of the expansion and therefore on the adiabatic losses, seems to be 
a "twin-corkscrew" or a simply S-shaped pattern (Hjellming & Johnston 1988; 
Crocker et al 2002). The last one seems to be the case of LS I +61° 303. Can we 
trace any ballistic motion of any jet component ? 

We have splitted the data base of each MERLIN run into two subsets. The 
first map (Fig. l37l-a^ represents the first four hours of the first run. It shows an 
ejection "A" already quite displaced from the core. In Figl37l-b. there is present 
a new ejection "B" at another PA. The combination of these two maps together 
produce as a consequence a "bent" jet, that is the southern jet of Fig.l35h. The two 
counter-jets for "A" and "B" are (Figs. l37l a and l37l-b) too weak to be detected, 
and they become visible only in the more sensitive image of Fig. l35k . where 
all the jet and counter jets for A and B form together the "S-shaped" jet. The 
"B" component is still detectable after 9 hours in the third image CFigl37l-c'). Its 
motion is ballistic: the PA is still the same (almost 90°) as in Figl37l-b. A new 
component "C" is present at another PA (Fig. [37j-d) 6 hours later, little rotation 
of the PA is compatible with A PA( B _c)/3 of the previous image. The Northern 
elongation in the higher sensitivity map of Fig.l35b therefore is the result of a set 
of little ejections of a rotating stream. The noise level in this image is still lower 
than that in Fig. l35k . nevertheless the counter-jet is not visible. This implies a 
decreased 9 due to precession. In the case of the MERLIN image of April 22 we 
derive [3cos0 = 0.12, which for (3 = 0.6 leads to an ejection angle of = 78°. 
This is an average of the ejection angles 0a and 0b of features A and B in Figs. 
l37k andl37b. A direct estimate of these angles is prevented because of the lack of 
the receding jets. Using the r.m.s. noise we derive 0a < 90°, 0b < 80° and for the 
C ejection in FigEZt, C < 68°. 

Therefore, the angle between the jet and the line of sight, 0, has decreased by 
more than 10° in 24 hours. It is this much narrower alignment of the jet with the 
line of sight, that causes the counter-jet to get further Doppler de-boosted with 
respect to the first image and lets it disappear below the sensitivity limit of the 
image. 
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Figure 36: : Morphology of a precessing and expanding jet. A precessing 
jet describes a twin-corkscrew pattern that projected in the sky plane resembles 
a quasi-sinusoidal path. Hjellming and Johnston (1988) have shown that only 
for slow expansion (Right-Bottom) it is possible to observe the jet in its twin- 
corkscrew pattern, whereas for free expansion (i.e. strong adiabatic losses) the 
morphology reduces to that on the Right-Top. The morphology in the center is 
that reproducing the observed structure (Bottom-left) and corresponds to a slow 
expansion switched to free expansion at a proper distance from the core. 



54 





100 




100 




100 


^ 






















D 









a) April 22 17h-22h 


-100 


b) April 22 22h- April 23 02h 


-100 


c) April 23 16h-22h 


-100 


d) April 23 22h- April 24 04h 



100 -100 100 -100 100 -100 100 -100 



Figure 37: MERLIN self-calibrated images of LS I +61°303 at 5 GHz using uni- 
form weights, obtained on 2001 April 22 and April 23. The data set of each epoch 
has been split into two blocks. A convolving beam of 40 mas has been used in all 
images for better display. The first contour represents the 3a level in all images ex- 
cept for c), where we start from the 2o level to display the faint B component. The 
rms noises are o=0.13 mJy beam -1 , o=0.20 mJy beam , a=0.13 mJy beam , 
and 0=0.15 mJy beam -1 , respectively. The PA of the ejections is indicated by a 
bar (see text). 

8 CONCLUSIONS 

The conclusions of this review of the astronomical methods used for the inves- 
tigation of Microquasars, with an examplary view on the source LS I +61°303, 
are: 

1. It is still an open issue, whether the compact object in this system is a neu- 
tron star or a black hole. In fact, taking into account the uncertainty in 
inclination, mass of the companion and the mass function, the existence of 
a black hole cannot be ruled out. 

2. The observational results from X-rays for LS I +61°303 are consistent with 
transitions between X-ray spectral states typical for a variable accretion 
disk. These transitions are properly related to the onset of strong radio emis- 
sion, as expected if the jet is "fed" by the disk. Quasi-periodic oscillations at 
soft X-rays and radio wavelengths are present in a strong analogy with those 
observed in GRS 1915+105. They occur at the onset and decay of large ra- 
dio outbursts. If confirmed, this fact might indicate that at the beginning the 
matter is ejected from the disk in the form of discrete condensations (i.e. 
blob-like), then follows a steady state where the matter-supply occurs at a 
higher/continuous rate (i.e. a continuous jet) and finally the ejection again 
ends up in a blob-like form. (Hujerat & Blandford 2003). 

3. At a scale of hundreds of AU the radio jet quite strongly changes its mor- 
phology in short intervals (< 24 hours), evolving from an initial double- 
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sided jet into an one-sided jet. This variation corresponds to a reduction of 
more than 10° in the angle between the jet and the line of sight. This new 
alignment severely Doppler de-boosts the counter-jet. Further observational 
evidence for a precessing jet is recognizable even at AU scales. 

4. The same population of relativistic electrons emitting radio-synchrotron ra- 
diation upscatters - by inverse Compton processes - ultraviolet stellar pho- 
tons and produces Gamma-ray emission. Ejections near the periastron pas- 
sage produce Gamma-ray flares but no radio flares, implying severe Comp- 
ton losses. 

We conclude that, because precession and variable Doppler boosting are the causes 
of the rapid changes in the radio-morphology, precession and variable Doppler 
boosting are likely to produce Gamma-ray variations at short time scales. The 
amplification due to the Doppler factor for Compton scattering of stellar photons 
by the relativistic electrons of the jet is S 3 ~ 2a (where a < 0), and therefore is 
higher than that for synchrotron emission, i.e. 8 2 ~ a (Georganopoulos et al. 2001; 
Kaufman Bernado et al. 2002). LS I +61°303 becomes therefore an ideal labo- 
ratory to test the recently proposed model for Microblazars with INTEGRAL and 
MERLIN observations now and by AGILE and GLAST in the future. 
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9 SUMMARY 



Because of their accretion disk super-massive black holes, with 10 6 -10 9 solar 
masses, in the heart of galaxies are the cause for the most energetic sources of 
emission in our Universe. The centers of such galaxies are called Active Galactic 
Nuclei (AGN). Some AGN, like the quasars, produce "jets" of subatomic parti- 
cles with speeds approaching that of light. A microquasar - as its name suggests 
- is a miniature version of a quasar: A disc of a few thousand kilometers radius 
surrounds a black hole of a few solar masses and two relativistic jets are propelled 
out of the disk by the same process occurring in a Quasar. The Microquasars 
therefore can serve as a convenient "laboratory" for studying the physics of jets. 
The Microquasars are objects very much closer to us than Quasars and the study 
of the evolution of relativistic jets can be done in a few days only, whereas on 
the contrary for far distant Quasars observations of many years apart are neces- 
sary to obtain appreciable proper motions of the radio jets. Moreover, concern- 
ing the intrinsic variability of Microquasars, these "small" objects change more 
quickly than Quasars: Considering as a characteristic time scale for variations 
t ~ Rschwarzschild I c Mass, phenomena of timescales of minutes connected with 
a Microquasar of 10 solar masses would take years in a AGN of 10 7 solar masses. 
Such an enormous difference is the main reason why Microquasars got such a 
great interest and growth in Astrophysics in the last decade. 

The Microquasars belong to the class of the X-ray binaries, where a compact 
object (black hole or neutron star) accretes from a normal companion star. Such 
systems are well known since the 1960s. The X-ray emission originates from the 
very hot accretion disk surrounding the compact object. However, it took a long 
time to discover that some of these systems also have relativistic radio jets like 
Quasars. For several years, after its discovery in 1979, SS 433 with its spectacular 
jets was thought to be a unique exotic case, a mere curiosity in our galaxy. Since 
the beginning of the 1990, after the discovery of other possible candidates of the 
same nature, several groups (including the author of this review) have begun a 
systematic research on X-ray binaries with radio jets. 

Here I review the astronomical methods used from Gamma-rays over X-rays 
and optical to radio wavelengths for the investigation of these objects. The de- 
scription of the methods is accompanied by directly applying them to the system 
LS I +61°303 , one of the most enigmatic objects in our galaxy, because it is as- 
sociated with a variable high-energetic Gamma-ray emission of unknown origin. 

The nature of the accretor - a neutron star or a black hole - is determined 
by optical measurements of the Doppler shift of spectral lines of the normal star 
orbiting around the invisible companion. Observations at X-rays probably are the 
most spectacular ones, in respect to the progress in the knowledge of the accretion 
disk and the disk-jet connection. Fitting the X-ray spectra information of the size 
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of the last stable orbit around the compact object can be derived and ejections of 
matter into relativistic jets can be related to variations of the disk. The results 
from X-ray observations for LS I +61°303 are consistent with transitions between 
spectral states typical for a variable accretion disk. These transitions are properly 
related to the onset of strong radio emission as expected for a jet "fed" by the disk. 

Onset and decay of some large radio outbursts are modulated with quasi- 
periodic oscillations that correspond to repetitive ejections of discrete conden- 
sations (i.e. blob-like). Continuous ejections have a flat radio spectrum and 
Low/Hard X-ray state. 

Speed and morphology of the ejections at high resolution are studied with ra- 
dio interferometric techniques. The results of more than 10 years of VLBI/EVN 
and MERLIN observations of LS I +61°303 are presented together with our dis- 
covery of the relativistic jet and its precession. Successive ejections are in ballistic 
motion; because of precession their projected path on the sky plane draw a bend- 
ing jet. 

The radio bursts occur around apastron passage, where the low velocity of the 
accretor enables it to capture more material of the wind from the companion star. 
However, no bursts are observed at periastron passage, where accretion theory 
predicts another super- accretion event. There the accretor is completely embed- 
ded in the densest part of the wind. I found that this great open question about 
LS I +61°303 and the other enigma about the association of LS I +61°303 with 
a variable Gamma-ray source are indeed not two separate questions, but on the 
contrary one is the answer to the other: We do not see a radio outburst at pe- 
riastron passage, because we see a Gamma-ray outburst. With LS 5039 for the 
first time we identified a Microquasar with an high-energy (E> 100 Mev) source. 
This fact opens the perspective that others of the more than one hundred still 
unidentified EGRET sources could belong to a new class of objects: Gamma-ray 
Microquasars. 

In this review I show that the variable Gamma-ray emission of LS I +61°303 is 
related to the orbit of the system, with peaks clustering where the companion star 
- a strong emitter of ultraviolet photons - is closest (at periastron). The suggested 
most probable explanation is that the ejected relativistic electrons are not able to 
emit synchrotron radiation at radio wavelengths, because at periastron passage 
they are embedded in such a strong UV-field of radiation that they loose com- 
pletely their energy by inverse Compton process. During the second accretion 
peak, the compact object is much farther away from the companion star and in- 
verse Compton losses are lower: The electrons can propagate out of the orbital 
plane and radio outbursts are observed. 
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10 ZUSAMMENFASSUNG 



Supermassive Schwarze Locher mit 10 6 -10 9 Sonnenmassen im Zentrum von Galax- 
ien sind wegen ihrer Akkretionscheibe der Grund fur die starksten Strahlungsquellen 
in unserem Universum. Die Zentren solcher Galaxien werden Aktive Galaktische 
Nuklei (AGN) genannt. Einige AGN wie z.B. Quasare produzieren "Jets" von 
subatomaren Teilchen mit Geschwindigkeiten bis nahe an die Lichtgeschwindigkeit. 
Ein Mikroquasar ist, wie der Name schon sagt, die Miniatur eines Quasars: Eine 
Scheibe von einigen tausend Kilometern umgibt ein Schwarzes Loch von eini- 
gen Sonnenmassen und zwei relativistische Jets werden durch denselben Prozess 
wie bei Quasaren aus der Scheibe herausgeschleudert. Ein Mikroquasar kann 
deshalb als ein brauchbares "Labor" zum Studium der Physik solcher Jets di- 
enen. Die Mikroquasare sind Objekte, die wesentlich naher zu uns liegen als 
Quasare und die Untersuchung der Evolution von relativistischen Jets kann in nur 
ein paar Tagen geschehen,wohingegen man fiir weit entfernte Quasare langjahrige 
Beobachtungen benotigt, um ausreichende Eigenbewegungen zu erhalten. Weiter- 
hin, wenn man die intrinsische Variability von Mikroquasaren betrachtet, andern 
sich diese "kleinen" Objekte schneller als Quasare: Nimmt man als charakteris- 
tische Zeitskala fiir Variationen x ~ Rschwarzschild l coc Mass, nehmen Phanomene 
mit einer Zeitskala von Minuten bei Mikroquasaren von 10 Sonnenmassen eine 
Zeit von Jahren bei einem Quasar von 10 7 Sonnenmassen in Anspruch, wenn man 
es mit der Masse des Akkretors skaliert.Diese enorme Differenz ist der Haupt- 
grund, weshalb Mikroquasare in der letzten Dekade solch ein grosses Interesse 
und Wachstum in der Astrophysik auf sich gezogen haben. 

Die Mikroquasare gehoren zur Klasse der Rontgen-Doppelsterne, wo ein kom- 
paktes Objekt ( Schwarzes Loch oder Neutronenstern ) von einem normalen Be- 
gleitstern einen Massenzuwachs erfahrt. Solche Systeme sind seit den 1960er- 
Jahren gut bekannt.Die Rontgenstrahlung stammt von der sehr heissen Akkre- 
tionsscheibe, die das kompakte Objekt umgibt. Allerdings dauerte es eine lange 
Zeit, bis man entdeckte, dass einige dieser Systeme ebenso wie Quasare relativis- 
tische Radio-Jets aussenden. Fiir viele Jahre nach seiner Entdeckung 1979 gait 
SS 433 mit seinen spektakularen Jets als ein einzelner exotischer Fall, eine einzi- 
gartige Kuriositat in unserer Galaxie. Anfag der 1990 nach der Entdeckung von 
anderen moglischen Kandidaten derselben Art begannen einige Gruppen ( ein- 
schliesslich des Autors dieses Reviews) mit einer systematischen Forschungsar- 
beit an Rontgen-Doppelsternen mit Radio- Jets. 

Ich gebe hier einen Uberblick iiber die astronomischen Methoden, die im 
Bereich von Gamma-Strahlung iiber Rontgen-Strahlung und optischen bis hin zu 
Radio Wellenlangen zur Untersuchung dieser Objekte angewandt werden. Die 
Beschreibung dieser Methoden wird unmittelbar begleitet durch die Anwendung 
der Methoden auf das System LS I +61°303 , einem der ratselhaftesten Objekte in 
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unserer Galaxie,weil es miteiner variablen hochenergetischen Gamma- Strahlungsquelle 
unbekannten Ursprungs verbunden ist. 

Die Natur des Akkretors - ein Neutronenstern oder eine Schwarzes Loch - wird 
durch optische Messungen anhand der Dopplerverschiebung von Spektrallinien 
des normalen Sterns, der sich um seinen unsichtbaren Begleiter bewegt, ermittelt. 
Beobachtungen im Rontgenbereich sind wahrscheinlich die spektakularsten im 
Hinblick auf den Fortschritt beziiglich der Kenntnisse iiber die Akkretionsscheibe 
und die Scheiben-Jet Verkniipfung. Indem man die Rontgenspektren untersucht, 
erhalt man Informationen iiber die letzte stabile Bahn um das kompakte Objekt 
und ebenso kann man die Variation dieser Grosse mit dem Auswurf von Materie 
in die relativistischen Jets correlieren. Die Ergebnisse von Rontgenbeobachtugen 
von LS I +61°303 sind konsistent mit Ubergangen zwischen spektralen Zustanden 
im Rontgenbereich typisch fur eine veranderliche Akkretionsscheibe. Diese Ubergange 
sind passend verbunden mit einem Anstieg von starker Radio strahlung, wie man 
es fur einen Jet, der von der Scheibe "gespeist" wird, erwartet. 

Anstieg und Abfall der starken Ausbriiche sind moduliert mit quasi-periodischen 
Oszillationen, die wiederholten Auswiirfen von diskreten Kondensationen entsprechen. 
Ununterbrochene Auswiirfen haben Radiospektrum flach und Rontgenzustand "Low/Hard". 

Die Geschwindigkeit und die Morphologie der Auswiirfe wereden mithilfe 
der Radiointerferometrie-Technik untersucht. Die Ergebnisse iiber mehr als zehn 
Jahre von Beobachtungen mit VLBJ7EVN und MERLIN von LS I +61°303 wer- 
den hier zusammen mit unserer Entdeckung des relativistischen Jets und seiner 
Prazession dargestellt. Aufeinander folgende Auswiirfe folgen ballistischer Be- 
wegung; wegen der Prazession bildet ihr auf die Himmelsebene projezierter Weg 
einen gebogenen Jet. 

Der Radioausbruch geschieht um den Apoastron-Durchgang herum, wobei die 
geringe Geschwindigkeit es dem Akkretor erlaubt, mehr Material vom Wind des 
Begleitstern einzufangen. Allerdings werden keine Ausbriiche beim Periastron- 
Durchgang beobachtet, wo die Akkretions-Theorie ein weiteres super-akkretives 
Ereignis vorhersagt. Dort ist der Akkretor vollstandig vom dichtesten Teil des 
Windes umgeben. Ich fand heraus, dass diese grosse offene Frage iiber LS I +6 1 °303 und 
das andere Ratsel beziiglich der Verbindung von LS I +61° 303 mit der veranderlichen 
Gamma-Strahlungsquelle in Wirklichkeit keine zwei getrennten Fragen sind, son- 
dern im Gegenteil ist die eine die Antwort auf die andere: Wir beobachten keinen 
Radioausbruch beim Periastron- Durchgang, weil wir einen Gamma-Strahlungsausbruch 
sehen. Mit LS 5039 haben wir zum ersten Mai einen Mikroquasar mit einer hoch- 
energetischen (E> 100 Mev) Quelle identifiziert. Diese Tatsache eroffnet die 
Perspektive, dass andere der mehr als hundert noch nicht identifizierten EGRET 
Quellen zu einer neuen Objektklasse gehoren konnen: Gamma-Strahlungs-Mikroquasare. 

In diesem Review zeige ich, dass die veranderliche Gamma-Strahlung von 
LS I +61°303 mit der Umlaufbahn des stellaren Systems verbunden ist, mit sich 
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laufenden Spitzen dann, wenn der Begleitstern - ein starker Strahler von Ultraviolett- 
Photonen - am nachsten ist (beim Periastron-Durchgang). Die vorgeschlagene 
wahrscheinlichste Erklarung ist, dass die ausgeworfenen relativistischen Elektro- 
nen nicht in der Lage sind, Synchrotonstrahlung im Radiobereich auszusenden, 
weil sie beim Periastron-Durchgang in solch einem starken UV-Strahlungsfeld 
eingebettet sind, dass sie vollstandig ihre Energie wegen des inversen Compton 
Prozesses verlieren. Wahrend der zweiten Akkretions-Spitze ist das kompakte 
Objekt wesentlich weiter weg von seinem Begleitstern und so sind Compton- 
Verluste geringer: Die Elektronen konnen aus der Bahnebene heraustreten und 
man beobachtet Radioausbriiche. 
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Figure 38: : Hybrid mapping. The loop "A" is the standard hybrid mapping 
scheme. Massi and Aaron (1999) have shown as the simple use of loop "B" elim- 
inates the spurious structures created by standard hybrid mapping. 

11 APPENDIX: Theory of Very-Long-Baseline-Interferometry 
Data Analysis 

The lack of phase information had prevented VLBI from being a true imaging 
technique until Rogers and his co-workers (1974) applied a phase closure rela- 
tionship. The introduction of the closure phase concept marks the beginning of 
a new era in VLBI. Many authors developed methods, reviewed by Pearson and 
Readhead (1984), which explicitly or implicitly use this quantity. Massi (1989) 
showed how the methods explicitly using the closure phase can be unified in one 
equation. In an attempt to unify all methods together Massi & Comoretto (1990) 
found that all methods turn out to be particular cases of the method proposed by 
Schwab (1980) depending on a proper scheme of baseline weighing. 

Using Schwab's method, called self-calibration,a map of the radio source can 
be obtained by using an algorithm which includes fourier transform and CLEAN, 
following an iterative procedure first indicated by Readhead and Wilkinson (1978) 
and called Hybrid mapping (Fig. l38b . That this procedure can converge on wrong 
solutions has been pointed out by many authors in the past: Walker (1986) in- 
dicated the bias in the resultant data due to the use of a point source as starting 
model. Baath (1989) suggested the use of the original data set in each iteration of 
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self calibration. Linfield (1986) analysed the role of the (u-v) coverage and lack 
of intermediate spacing. Generally, the full procedure to avoid false features is not 
clear and only experience with imaging helps the user to avoid them. 

Mas si & Aaron (1999) demonstrated that the problem is connected with the 
non-linear nature of self calibration which leaves initial wrong assumptions frozen 
in the final solution. We demonstrated that the general precondition to avoid false 
structures in the map is that the errors (or more precisely their cube) of the model 
should be smaller than the observed closure phases. This condition, generally 
satisfied for a standard earth based array, is violated if one telescope of the array is 
very displaced from the others, as it is for an array including a telescope mounted 
on a satellite. In this case one should avoid the use of a point like model as starting 
model. Moreover, one should at each iteration of self calibration, adopt the model 
derived by CLEAN directly on the original data and not on the corrected data 
biased by previous wrong solutions (Fig. l38t . 

Using self calibration it is assumed that the baseline based errors are negligi- 
ble. In spite of the fact that these baseline errors are quite small their effect on the 
map's quality is rather serious. Tests were performed to determine at which level 
errors limit the obtainable dynamic range with the VLA (Perley 1986), with the 
VLBA (Briggs et al. 1994), with the Nobeyama Radioheliograph (Koshiishi et 
al., 1994). Massi and collaborators have performed such an analysis for the Euro- 
pean VLBI Network (EVN). The result was that the instrumental polarization (D 
terms) of the telescopes of the network had an average value of 9 percent arriving 
at some telescopes at values of 20 percent. For comparison the values of VLBA 
telescopes were below 2 percent. The instrumental polarization was therefore the 
main reson for the lower performance of the EVN damaging the dynamic range of 
the images up a factor of 7 (Massi et al. 1991; 1996; 1997; 1997b; 1997c; 1997d; 
1997e; 1998; Massi 1999b). 
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14 ABSTRACT 



The Astrophysics of microquasars - galactic miniatures of the far distant quasars 
- has become one of the most active fields of modern Astronomy in recent years. 
Here I review the astronomical methods used for the investigation of these objects, 
from Gamma-rays over X-rays and optical to radio wavelengths. The description 
of each astronomical method is always followed by an examplary application on 
the source LS 1+61 303, one of the most observed Be/X-ray binary systems be- 
cause of its periodical radio emission and strong, variable Gamma-ray emission. 
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